Arbitrary-order non-linear contribution to self-steepening 
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Abstract 

Based on the recently published generalized Miller formulae, we derive the spectral dependence 
of the contribution of arbitrary-order non-linear indices to the group-velocity index. We show that 
in the context of laser filamentation in gases all experimentally-accessible orders (up to the 9 th - 
order non-linear susceptibility hi a h an d x'' 11 ^ i n argon) have contributions of alternative signs 
and similar magnitudes. Moreover, we show both analytically and numerically that the dispersion 
term of the non-linear indices must be considered when computing the intensity-dependent group 
velocity. 

PACS numbers: 190.3270 Kerr effect; 190.7110 Uitrafast noniinear optics; 120.4530 Opticai constants 
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Non-linear optics jl| relies on the non-linear properties of the propagation medium, among 
which the successive orders of the non-linear susceptibility are essential parameters. How- 
ever, due to the difficulty to measure them experimentally, their knowledge is generally 
limited to the first non-zero order (x^ or X ■> depending on medium symmetry). Fur- 
thermore, the available laser sources drastically limit the wavelengths available for such 
measurements, so that reliable dispersion curves for higher-order susceptibilities cannot be 
deduced from the sparse experimental data available to date. 

The lack of data led to neglect these higher-order Kerr terms in most numerical simula- 



tions of e.g. both self-guided filaments in ultrashort intense laser pulses 



5] or the propa- 



gation of high-intensity pulses in hollow-core fibers jfj]. Similarly, the Kerr contribution to 
the group velocity is most generally limited to the third-order and treated as dispersionless 
in the lack of data about its dispersion Jjj, Q]. Recently, however the measurement of the 
higher-order refractive indices up to n 8 in No and 2 , and up to n w in argon followed 
by the generalization of the Miller formulae js| to any order of non-linearity [9], provided a 
new insight into the spectral dependence of the non-linear refractive index at high incident 
intensity. 

This allowed us to show that these terms cannot be ne glec ted and can even provide the 
dominant contribution stabilizing the self-guided filaments 10|. Furthermore, in argon-filled 
hollow-core fibers, these terms are necessary to obtain quantitative agreement of numerical 
simulations with experimental data III, [12] . But these works focused on the contribution 
of higher-order Kerr terms to phase velocity. The contribution of the spectral dispersion of 
higher-order indices to group velocity was not considered, although it can be expected to 
impact the propagation, and in particular the self-steepening term. 

Here, we derive an explicit expression for the contribution of any order of the non-linear 
refractive indices to the group-velocity index. We show that in filamentation, all non- 
linear orders of the group- velocity index can have similar orders of magnitude and must 
be considered. Furthermore, especially in the ultraviolet, the dispersion of the non-linear 
indices cannot be neglected when calculating self-steepening. 

At arbitrary intensity / and frequency u, the refractive index can be expressed as 

oo 

n(uj) = n (u) + n 2 {uj)I + n A (u)I 2 + ... = s ^n 2 j(u)P (1) 

j=0 
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In gases, where n_^- 1 < 1, the refractive index n(u) = 1 + Y^T=o X^ +1 K U} )I-' can be 



approximated by 



Q: 



n (co) « l + (2) 
Z W+i) x (^+i)( w ) ( 3 ) 

where x^ +1 ^ is the (2j + l)' ft -order non-linear susceptibility and the Z^ 2 ^ +l "> are frequency- 
independent factors. If a pulse can be described as a carrier wave modulated by an envelope 
with a sufficiently narrow spectrum to allow neglection of the envelope deformations over 
short distances, then a group-velocity index can be defined as 

nJuj) = — = n(tu) + u^- (4) 

Vg Q.UJ 

where v g is the group- velocity and c is the speed of light in vacuum. Defining = 1/2 
and considering Equations (JT])- (J3j) , n g rewrites: 

00 / /-/ (2.7+1) \ 

« i + ^ (2i+1) ( X {23+1) (w) + ^ ) (6) 

Identifying the terms for each power of the intensity, we obtain the contribution of each 
order of non-linearity to the group- velocity index. 



V 3 > 1, „ S , 23 ( W ) „ Z««> (>«> + (8) 

The first expression corresponds to the usual linear contribution, while the first non-linear 
contribution to the group- velocity index n 2j9 corresponds to the classical self-steepening term 

II HQ. 

Within the elastically bound electron modeljthe susceptibility of abitrary order 
is given up to any order by generalized Miller formulae |9[: 

me {(jJo z - u 2 + «w 7) 

vg>2, = — (-) fl Q (9, o7^r ( 10 ) 



where m and — e are the electron mass and charge, eo the permittivity of vacuum, N the 
density of dipoles in the propagation medium, 7 the width of the resonance at frequency ujq 
and describes the potential well where the electron oscillates; Q(oj) = to 2 — to 2 + icoo'j. 
Inserting these expressions into Equations (j7|) and (JSJ) yields: 

n 9i oM-l - M^)-l) ^_^. o7 (11) 
V, > l,n fl , 2 » « n^M ^ + ^+y + fcvy (12) 

This equation provides a general expression of each non-linear contributions to the group- 
velocity index n g . As a consequence, it allows to evaluate the impact of higher-order Kerr 
terms on the self-steepening in the context of laser filamentation in gases or the propagation 
of ultrashort pulses in hollow fibers. In the following, we consider the propagation of high- 
intensity pulses in transparent media, far from resonance. In this case, \ojq — oj\ ^> 7, so 
that Uq + u 2 > |o>q — u 2 \ ^> cj 7. The imaginary parts of Equations ffTTj) and f[T2l become 
negligible: 

2 2 

%oM-i « KM-i)^ (13) 

. . u 2 + (4j + 3)w 2 

WW) « n 2i 2 / j2 ( 14 ) 

Note that the negative values obtained for u; > cjq correspond to the well-known region of 



negative group- velocity index 15j. Figured] displays the spectral dependence from Equation 



( I14p . based on the recent experimental measurements of n 2 j at 800 nm 
to the whole visible spectrum by applying generalized Miller formula? 



16]. 
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10|, extrapolated 



and the dispersion 



data of Zhang et al. 

From Equations ( 1T31) and ( JT41) . we can estimate the ratio of the successive terms of the 
group- velocity index 

n g ,2l n 2 I 0Oq + 7co 2 



n gfi - 1 n - 1 u 2 + u 2 
n g , 2j+2 P +1 _ n 2j+2 I ou 2 + (4j + 7)u 2 



(15) 



vj > 1, g '^;, = r. : 2 ae) 



w 9 ,2i/ J ' n 2 j w 2 + (4j + 3)w 2 



The second factor of equations (fl5l) and ( fT6l) is of the order of 1. Therefore, the orders of 
magnitude of the ratio of successive terms are driven by the ratio of the non-linear indices 
n 2 j, multiplied by /. The values displayed in Figure []] imply that, for / < 10 14 W/cm 2 , 




FIG. 1: Spectral dependence of the non- linear group- velocity indices (a) n g< 2, (b) n 9i 4, (c) n 5i 6, 
and (d) n ffj g of O2, N2, air and Ar at atmospheric pressure 

n 2jP ^ n — 1. Self-steepening is therefore, as well-known Q], a second-order term in 
the non-linear Schrodinger equation (NLSE) describing the non-linear propagation of light 
in a non-linear transparent medium. Furthermore, all known terms in n-ijP have alternate 
signs and comparable orders of magnitude . The same therefore applies to the terms 
in n g> 2jP, which must all be taken into account when describing self-steepening e.g . in the 



context of filamentation, where the intensity is clamped around 5 x 10 13 W/cm 2 17|, [l8|, or 
of the propagation in hollow fibers. 

Equation (1141) provides an estimation of the error performed when neglecting the disper- 
sion term in the contribution of the higher-order indices to the group- velocity index. Figure 
[2] displays the relative error 1 — n 2 j/n gt 2j implied at atmospheric pressure when neglecting 
the dispersion terms of the Kerr contributions to the group velocity. The calculations are 
based on the same data as in Figure [TJ As is clear from Equation (JHJ), this error decreases 
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FIG. 2: Relative error induced when neglecting the dispersion of the Kerr terms in the group- 
velocity (a) 1 — n2/rig t 2, (b) 1 — n^/ng^, (c) 1 — n§/n g fi, and (d) 1 — n%/n g $ of O2, N2, air and Ar 
at 1 atm. 

for longer wavelengths where dispersion is smoother. At 800 nm, it amounts to ~20% 
and may be considered acceptable, although not negligible. However, at blue or ultraviolet 
wavelengths, the dispersion term dominates and must be considered in the equations. 

Numerical simulations of the propagation of a 35 fs pulse in a 1 m long hollow-core 
fiber filled with 1.4 bar argon, confirm this finding. As described in detail earlier 
model implements the NLSE including the higher-order Kerr terms. We compared the code 
output with and without the contribution of the higher-order indices to the group-velocity 
index up to the term in nio, i.e. the terms of Equation (|6]) for 1 < j < 5. As can be seen in 
Figure El the consideration of the full steepening term affects the spectrum by deforming the 
pulse envelope. It simultaneously red-shifts the central part of the spectrum, and blue-shifts 
its edges. Furthermore, as predicted by the analytic calculations, the contribution of the 
dispersion of the higher-order Kerr terms is larger in the UV and negligible in the infrared. 



12|, the 
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FIG. 3: Influence of the higher-order indices contribution to the group- velocity index on the prop- 
agation of a 35 fs, pulse with fixed 712(A)/ in a 1 m long hollow-core fiber filled with 1.4 bar argon. 
(a,d) 250 nm, 286 //J; (b,e) 400 nm, 360 /xj; (c,f) 800 nm, 400 fri 

These these terms must therefore be considered in numerical simulations, especially while 
investigating spectral broadening. 

In conclusion, based on the recent generalization of the Miller formulae, we have estimated 
the contribution of higher-order indices to the group- velocity index. These contributions de- 
fine the self-steepening term of the NLSE. They have alternate signs and comparable absolute 
values in intensity regimes typical of filament at ion. All non-linear terms must therefore be 
considered in the evaluation of the self-steepening of ultrashort intense laser pulses prop- 
agating in transparent Kerr media. Furthermore, we demonstrate both analytically and 
numerically that their spectral dispersion cannot be neglected either, especially at shorter 
wavelengths. 
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